Preliminary lattice QCD calculation of the a meson decay width 
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We present an exploratory lattice QCD calculation of the a meson decay width using the s-wave 
scattering phase shift for the isospin 7 = pion-pion (tttt) system. The finite size formula developed 
by Rummukainen and Gottlieb is employed to estimate the phase shift, which demonstrates a 
characteristic consistent with the presence of a resonance around a meson mass. The effective range 
formula is adopted to describe the phase-shift dependence, and we extract the effective cr — >■ tttt 
coupling constant Qa-ww = 2.69(44) GeV, which is consistent with the theoretical predictions. The 
decay width estimated from phase shift is about 236 ± 49 MeV. Our preliminary lattice simulations 
are carried out with MILC full QCD gauge configurations in the presence of 2 + 1 flavors of the 
"Asqtad" improved staggered dynamical sea quarks on a 16'^ x 48 lattice at ruTi/ma ~ 0.414 and 
lattice spacing a « 0.15 fm. 
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I. INTRODUCTION 

It is well-known that sigma meson is a resonance, which 
is a state with a considerable decay width under strong 
interactions, and decay via strong interaction. In 2010, 
Particle Data Group (PDG) [J lists the meson /o(600), 
which is usually called a meson [J^^ = 0++, / = 0), 
with an estimated mass (400 — 1200 MeV) and an esti- 
mated decay width (600 - 1000 MeV). Such a low-lying 
scalar meson has been known to play a significant role in 
nuclear and hadron physics. The existence of the a meson 
has been established by some refinement of experimental 
analyses [2, jQj and some phenomenological studies [l2| - 
[2l| . The Dalitz plot analysis of Fermilab experiment 
E791 yields its decay width about 324^^^ ± 21 MeV. 
Moreover, the a meson has been extensively studied with 
BES data [1, 01 , and the most recent analysis based on 
the events collected by BESII gives the pole position of 
to be cr (541 ± 39) + i{252 ± 42) MeV [|. 

Although the direct determination of the a resonance 
parameters from QCD is afflicted with the fundamen- 
tal difficulties since the calculation of resonance masses 
and decay widths is essentially a nonperturbative prob- 
lem. Several research groups have undertaken theoretical 
efforts to study the a meson and its resonance parame- 
ters p^ - [2l| . Originally the cr-meson was introduced to 
fit experimental data and its mass was chosen to repro- 
duce the experimental results. There is a wide variety 
for defining the mass and width of a particle. Some au- 
thors use the pole approach with the mass and width of 
resonance taken from the position of the pole of the T- 
matrix p^ . Another way to study the mass and width 
of resonances is through the use of the spectral func- 
tion [2^. However, all the experimental and theoretical 
analyses give a little bit different resonance mass and de- 
cay width for a meson. Therefore, there is not close to 
the consensus yet on its resonance parameters. 

The most feasible way to extract a resonance parame- 
ters nonpcrturbatively from first principles is the use of 
lattice QCD. The dominant a decay channel is to a pair of 



pions, which can be precisely handled on the lattice. Un- 
til now, no direct lattice QCD study about a resonance 
have been reported, possibly because the rectangular and 
vacuum diagram are extremely difficult to reliably calcu- 
late. Inspired by J. Nebreda and J. Pelaezet's theoretical 
studies on a resonance [2^ , o ur p revious works about the 
precise extraction of a mass [23| , our reliable extraction 
of the TTTT scattering length at the 7 = channel [2^ . 
and our preliminary lattice QCD calculation of k meson 
decay width [2H , [26| , in this work, we further explore its 
resonance parameters directly from lattice QCD. 



In the present study, we estimate the a decay width 
by calculating the s-wave scattering phase shift for the 
/ = TTTT system. We discuss the ground state of the 
1 = TTTT system with the total zero momentum in cen- 
ter of mass (CM) frame, and the non-zero momentum 
in the moving frame (MF). The calculations are car- 
ried out with MILC full QCD gauge configurations in 
the presence of 2 -I- 1 flavors of the Asqtad improved 
staggered dynamical sea quarks, generated by the MILC 
collaboration [23, Hsj. The meson masses were deter- 
mined from our previous spectrum analysis p3| which 
gave rriT^/mc ~ 0.414, and the lattice parameters were 
determined by the MILC collaboration, namely, the lat- 
tice extent L is about 2.5 im and the lattice space inverse 
l/a = 1.358 GeV [13, [11]. The finite size formula [H- 
|33| is employed to the tttt system in the moving frame, 
and we utilize the Rummukainen- Gottlieb formula (3^ . 
which is the extension of the Liischer formula to the mov- 
ing frame, to estimate the s-wave phase shift. This cal- 
culation is carried out at two energies which enable us to 
investigate the existence of the a resonance. 

This paper is organized as follows. In Sec.|ll]we discuss 
our calculation method. In Sec. IIIIl we show the simu- 
lation parameters and our concrete lattice calculations. 
We present our lattice simulation results in Sec. II VI Our 
conclusions and outlooks are given in Sec. [V] 
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II. METHODS 



A. Scattering phase 



In the TTTT system, the relativistic Brcit-Wigner formula 
(RBWF) for the elastic s-wave scattering amplitude ao 
in the resonance region with a center-of-mass energy 
and a decay width Tji reads [l| 



ao 



il + iy/sT]i{s) 



(1) 



where s = E'^,^ is the invariant mass of the tttt system, 
EcM is the center-of-mass energy, and ao is related to the 
s-wave scattering phase So through ao = (e^*''" — l)/2i. 
The RBWF can be conveniently expressed as 



tan 5q 



(2) 



The a resonance has quantum numbers I{J^) = 0(0+) 
and decays into two pions in the s-wave. The decay width 
r7}(s) can be written in terms of the effective a — ^ tttt 
couphng constant gaTr-ir [13, 



p 

TT s ' 



(3) 



where p = ■\/s/4 — is the center-of-mass momentum 
of the pion. With a combination of Eqs. ([2]) and (jSj, a 
description of the scattering phase shift in the s-wave as 
a function of the ^/s is provided by so-called the effective 
range formula (ERF) in the elastic region, namely. 



9(T7V7r 



8tt \/s{ml - s) 



(4) 



This equation allows us a linear fit or solving for two 
unknown parameters: coupling constant g^T^Tr and the 
resonance position m^. 

It should be stressed that we will use the ERF to solve 
the scattering phase shifts which are measured by our 
lattice simulations, where we even used the pion mass 
which arc larger than its physical ones. 

The a decay width Fo- can then be calculated through. 



rfl(s) 



9o 



Pa 



8tt 



Pa 



Therefore, Eqs. ^ and (O allow us to extract mc and 
T„ through the dependence of the s-wave tttt scattering 
phase shift Jo on the invariant mass s. 



B. Finite- volume method 

A calculation of the phase shift from lattice QCD is 
practicable by using the finite-size formula proposed by 



Liischer |29l - l33l |. In this method, the phase shift is ex- 
tracted from the energy eigenvalues of the tttt system 
enclosed in a cubic box with spatial extension L. In this 
work, we focus on the a meson decay into two pions in 
the s-wave, and only consider the tttt system with the 
isospin representation of {I,Iz) = (0,0). 



1. Center of mass frame 

In the center-of-mass frame, when the a meson is at 
rest, the possible energy eigenvalues of the noninteracting 
TTTT system are given by 



E^2^ml+p\ (6) 

where p = |p|, p = (27r/L)n, and n E Z,^. In a typical 
lattice simulation, this energy for n 7^ is much larger 
than the sigma mass rria-. For example, on our full QCD 
configurations, the lowest energy for n ^ calculated 
from 7Ti7r and rua determined in our previous study [23| . 
is E = 1.56 X m^r, which is obviously away from the 
resonance, and thus it is not eligible to study the cr meson 
decay. Thus, we only consider the dominant low energy 
states, namely, n = case, where pions and a meson are 
at zero momentum. The energies of these states are 



El 
E2 



2mT, 

nirr 



for TTTT system, 
for (7 meson . 



(7) 



On our full QCD configurations the invariant mass of two 
free pions takes -^s = 0.828 x rUa, which is still a little 
bit away from the resonance mass mg., but closer to rric 
than those with n 7^ 0. 

In the interacting case, the total momentum of the tttt 
system is also zero. The energy eigenvalues are shifted 
by the hadronic interaction from En to En {n = 1,2), 
and the energy eigenvalue for tttt system is given by. 



E = 2^r 



k = 



27r 



(8) 



(5) ^ 



where q is no longer required to be a integer. In this 
work, we restrict on the energy cigenstates with energies 
E in the elastic region 2m,r < E < 4m7r with tttt system 
having the same quantum numbers as a meson. In the 
center-of-mass frame these energy eigenstates transform 
under the cubic group Oh in the irreducible representa- 
tion F = . The finite-size formula rela ting the energy 
En to the scattering phase shift So is [29l - l33l | 

tanSoik) = — , 
2oo(l;'7 ) 

where the zeta function denoted by 

1 x-^ 1 



(9) 



2oois;q ) 



E 



(10) 



The expansion of the zeta function in Eq. (jlOp is con- 
vergent only when Re s > 3/2, but it can be analytically 
continued to s = 1. Zoois; q^) can be efficiently evaluated 
by the method described in Appendix A of Ref. [s^ . 
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2. Moving frame 

In order to realize the energy of the tttt system is more 
close to the sigma mass nia--, we consider a system hav- 
ing a non-zero total momentum, namely, the moving 
frame [33 |. Using a moving frame with total momen- 
tum P = (27r/L)d, d G Z'^, the energy eigenvalues for 
non-interacting tttt system are 



Emf — 



Pi 



(11) 



where pi = |pi|, and pi, P2 denote the three-momenta 
of the pions, which satisfy the relations. 



—Ill , rij e Z'^ 

1j 



Pi + P2- 



(12) 



In the center-of-mass frame, the value of the Ecm is 
directly given by the discrete energy eigenvalue extracted 
from the large time behavior of the corresponding corre- 
lation function. The rclativistic 4-momentum squared 
is invariant, and Ecm is related to Emf in the moving 
frame through the Lorentz transformation 



E. 



CM 



E 



MF 



(13) 



In the moving frame, the center-of-mass is moving with 
a velocity of v = P/Emf- Using the standard Lorentz 
transformation with a boost factor 7 = l/Vl — v^, the 
Ecm can be obtained through 



Ecm = I^'Emf - 2^ml+p*^, 



(14) 



where, in the center-of-mass frame, the total center-of- 
mass momentum vanishes, namely. 



P 



= -P2 = ^7 ^(Pi -P2), (15) 



here and later we denote the center-of-mass momenta 
with an asterisk (*), the boost factor acts in the direction 
of V, and we use the shorthand notation 



7 ^P 



7 P|| 



where p|| and p_L are the components of p parallel and 
perpendicular to v, namely, p|| = (p-v)v/|vp, and pj_ — 
p — P||. From inspecting Eqs. (fT2)) and (fTS)) . wc can note 
that the p* arc quantized to the values 



. 27r 



rePd 



where the set is 



Pd= {r r = 7" 



n e Z-^ 



(17) 



(18) 



In this work, we only consider the dominant low en- 
ergy states in the moving frame, which are one pion 
with zero momentum, and one pion with the momentum 



p = (27r/L)e3 (d = 63) and a meson with the momentum 
P = p. The energy eigenvalues of these states are 



El = + y/m^ -I- 
E2 - ^rnl + P2 



for TTTT system, 

for a meson, (19) 



where P — |P|. On our full QCD configurations the 
invariant mass of the tttt system takes y/s = 0.994 x rricr, 
which is much closer to a mass nia- than that with total 
zero momentum in the center-of-mass frame. Therefore, 
we will only consider this case. The invariant mass ^/s 



was evaluated from ^/s = y/E^ 



MF 



P2. 



In the interacting case, the Ecm is given by 



k = ^q. 



Ecm = 2^ ml + 



(20) 



where q is no longer required to be a integer. It is ex- 
actly this energy shift between the non-interacting situa- 
tion and the interacting case, namely, Ecm — Ecm, that 
we can extract the tttt scattering phase. In this work, 
we only consider one moving frame, namely, d = 63, 
where the energy eigenstates transform under the tetrag- 
onal group -D4/1. The irreducible representation A2 is 
relevant for the tttt scattering states in the infinite vol- 
ume. Wc calculate the energies associated with A2 sec- 
tor. The hadron interaction shifts the energy from En to 
En {n ~ 1,2), and the energies En are related to the tttt 
scattering phase shift So in the infinite volume through 
the Rummukaincn- Gottlieb formula |34| . namely. 



tan(5o(fc) = 



77r^/2g 



(21) 



where we assume that the phase shifts di with I = 
2,4,6,... are negligible in the energy range of interest, 
and the zeta function is formally defined by 



1 



rePd 



(22) 



and the set Pd is denoted in Eq. ([T5)) . The zeta function 
converges when Re 2s > 1+3, and can be analytically con- 
tinued to whole complex plane. The k is the momentum 



(16) defined from the invariant mass -^/s as ^/s = 2y/k? 



The calculation method of Z^q{1; q^) is discussed in Ap- 
pendix A of Ref. nil . We also discussed it in Ref. [23 , 
where we extend it to more general case. Using Eq. <\22\ , 
we can obtain the phase shift from the energy eigenvalue 
calculated in the lattice simulations. 



C. Correlation matrix 

In the center-of-mass frame, the value of the Ecm is 
directly obtained by the discrete energy eigenvalue ex- 
tracted from the correlation function at large time t. 
Ecm is related to through the Lorentz transfor- 

mation 



Ecm 



^mf 



(23) 
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To calculate two energy eigenstates, i.e., En [n = 1,2), 
we build a matrix of the time correlation function, 



C{t) 



(o|oLWa.(o)|o) (o|oLWa(o)|o) 

(0|Ot(i)o^^(0)|0) (0|Ot(i)o^(o)|o) 



(24) 

where [t) is an interpolating operator for the a meson 
with the definite momentum, and 0^7^(0 interpolat- 
ing operator for the tttt system with specified momentum. 
The interpolating operators Oa{t) and O^TrCi) employed 
here is exactly the same as in our previous studies in 
Refs. [i^, m, li^l- The notations adopted here are also 
the same, but to make this paper self-contained, all the 
necessary definitions will be also presented below. 



1. TTTT sector 

In this subsection, we briefly review the formulas of 
the scattering length used in the present work from two- 
particle energy in a finite box, with emphasis on the for- 
mula for the / = s-wave tttt system. Also we discuss 
the detailed procedure for extracting the energies of the 
TTTT system. Here we follow the original derivations and 
notations in Refs. [sB-Hll- 

Let us consider the tttt system of two Nambu- 
Goldstone pious in the Asqtad-improved staggered dy- 
namical fermion formalism. Using operators Ot,{xi), 
OTr{x2), OTr(x^), OTr{x4) for pious at points xi, X2, x^ 
and X4, respectively, with the pion interpolating field op- 
erators defined by 



n-it) 
-'it) 



-^d(x,i)75it(x,i) 

X 

^it(x,i)75d(x,t) 

X 

*)'^5w(x, t) - u(x, t)75d(x, t)]. 



We then represent the tttt four-point functions as 

C7r,r(a;4,a;3,x2,a;i) = {OTr{xA)0-„{x3)OX{x2)0'l{xi)) , 

where (■ ■ ■ ) represents the expectation value of the path 
integral, which we evaluate using the lattice QCD simu- 
lations. 

In the current study, we only consider the dominant 
low energy states in the moving frame, which are the 
one pion with zero momentum, and one pion with the 
momentum p = (27r/L)e3 (namely, d = 63) and the 
total momentum P = p. The operator which creates a 
single pion with non-zero three momentum k from the 
vacuum is obtained by the Fourier transform: 

0,(k,i) =^a(x,t)e'''-^. (25) 

X 

After summing over the spatial coordinates xi, X2, X3 



and X4, we obtain the tttt four-point function in the 
mentum p state. 



mo- 



C'ir7r(P, ^4, ^3: ^2, il) — 

E 

Xi,X2,X3,X4 



e'P-("^-"^)a.(a;4,a;3,X2,a:i),(26) 



where xi = (xi,ti), X2 = (x2,t2), 2:3 = (x3,t3), and 
2:4 = (x4, ^4), and t stands for the time difference, namely. 

To avoid the complicated Fierz rearrangement of quark 
lines, we used creation operators at time slices that are 
different by one lattice time spacing as it is suggested in 
Ref. [11], namely, we select ti = 0,^2 = 1,^3 = and 
U^t + l. 

In the TTTT system, there are two isospin eigenstates, 
namely, / = and / = 2, in this study, we only consider 
/ = channel, we construct the tttt operators for this 
isospin eigenchannel as 



O 



^-"(p,t) = -^{^-(t)^+(p,i + l) + 7r+(i)^-(p,t + l) 
-^°(V(P,i + l)}, (27) 



where the p is the total momentum of the tttt system 
or the momentum of pion meson. This tttt operator has 
(/,/,) = (0,0). 

If we assume that the u and d quarks have same mass, 
only four diagrams contribute to the tttt scattering am- 
plitudes j36l - l38l |. The quark line diagrams contributing 
to TTTT four-point function are displayed in Fig.[Tl labeling 
them as direct(£'), crossed (C), rectangular (i?), and vac- 
uum {V) diagrams, respectively. The direct and crossed 
diagrams can be easily evaluated by constructing the cor- 
responding four-point amplitudes for arbitrary values of 

and ^4 using only two wall sources placed at the fixed 
time slices tx and t2- However, the rectangular (i?) and 
vacuum diagrams {V), require another quark propaga- 
tors connecting the time slices ^3 and t4. 




(a) D (b) C 



FIG. 1: Diagrams contributing to tttt four-point functions. 
Short bars stand for wall sources. Open circles are sinks for 
local pion operators. 

Encouraged by the exploratory works of the tttt scatter- 
ing at the / = channel [s^, , and our reliable extrac- 
tion of the TTTT scattering length [2J] at the / = channel 
with zero momentum, in the same way, we solve this 
problem by evaluating T quark propagators on a L'^ x T 
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lattice in the same way, each propagator, which corre- 
sponds to a wall source at the time slice t = 0, - ■ ■ ,T —I, 
are denoted by 

Y,Dn'.n"Gt{n") = Y,Sn',i^,t), 0<i<T-l, (28) 

n" X 

where D is the quark matrix for the staggered Kogut- 
Susskind quark action. The combination of Gt{n) that 



we apply for tttt four-point functions is shown in Fig. [T] 
For the non-zero momentum, we used an up quark source 
with 1, and an anti-up quark source with e^^'^ (except 
for D, where we use 1) on each site for two pion creation 
operator, respectively. D, C , R, and V , are schematically 
shown in Fig.[Tl and we can express them in terms of the 
quark propagators G, namely. 



G^Trip, t4, <3, <2, ^l) 
^L-(P; ^4,^3,^2,^1) 



I^^EE (Tr[G,\(x3,t3)Gt,(x3,i3)Gt(x4,t4)Gt,(x4,i4)] 

X3 X4 

I^^EE e*''"' (Tr[Gl^(x3,t3)Gt,(x3,i3)Gljx4,t4)Gt,(x4,i4)] 



X2 X3 

i^^EE' 



,JPX2 



Tr[Gj^ (x2, t2)Gt, {^2M)Gl (X3, h)Gt, (X3, ^3)] 



.ip-|,X2-X3) 



Tr[G,\(x2,t2)Gt,(x2,t2)G^X3,t3)Gu(x3,t3)]) , (29) 

I 



where daggers mean the conjugation by the even-odd par- 
ity (—1)" for the staggered Kogut-Susskind quark action, 
and Tr stands for the trace over color index. The her- 
miticity properties of propagator G are used to eliminate 
factors of 7^. We know that when we calculate the dis- 
connected part of the sigma correlator for the cases with 
non-zero momenta, the subtraction of the vacuum expec- 
tation value is not needed in Refs. [l^, [s^H^I. Similarly, 
the vacuum diagram here is not accompanied by a vac- 
uum subtraction for the cases with non-zero momenta. 

As it is discussed in Refs. [13, H^, the rectangular and 
vacuum diagrams create gauge-variant noise. One can 
reduce the noise by fixing gauge configurations to some 
gauge ( e.g.. Coulomb gauge), and select a special wall 
source to emit only the Nambu-Goldstone pion [4l|, how- 
ever, the gauge non-invariant states may contaminate the 
TTTT four-point function. Alternatively, we perform the 
gauge field average without gauge fixing since the gauge 
dependent fluctuations should neatly cancel out among 
the lattice configurations. Besides these cancelations, the 
summation of the gauge-variant terms over the spatial 
sites of wall source further suppresses the gauge- variant 
noise. In this work we found that it works well. 

All four diagrams in Fig. [1] are needed to study tttt 
scattering in the / = channel. As it is investigated in 
Ref. (37| . in the isospin limit, the tttt correlation function 
for J = channel can be expressed as the combinations 
of four diagrams, namely, 

C^^{t) = (C(t)IC(O)) =D+^C-iR+^V , (30) 

where the operator O^^t^ denoted in Eq. (I?7l) creates a tttt 
state with total isospin 0. 

In the usual manner, pion mass can be evaluated 



through one-pole fitted functions 

C^{t) = cosh [-m^{t - T/2)] + 



(31) 



In our concrete calculation we also evaluate the ratios 



i?^(0 



G^,(p,0,l,^,t + 1) 
G.(O,0,i)G.(p,l,t + l)' 



X= D,C,R, andV, 
(32) 

where Gt(0, 0, t) and G^(p, 1, t+1) are the pion two-point 
functions with the momentum and p, respectively. 

We should keep in mind that, for the staggered Kogut- 
Susskind quark action, there are further complications 
in itself stemming from the non-degeneracy of pions 
in the Goldstone and other channels at a finite lat- 
tice spacing. Briefly speaking, the contributions of non- 
Nambu-Goldstone pions in the intermediate states is ex- 
ponentially suppressed for large times due to their heav- 
ier masses compared to these of the Nambu-Goldstone 
pion j36l - l38| . Therefore, we consider such tttt interpola- 
tor does not couple significantly to other tttt tastes, and 
neglect this systematic errors for the tttt sector. 



2. a sector 



In our previous works [2j, |3^, |40l , we make a detailed 
procedure to measure (0|iT'''(t)(T(0)|0). To simulate the 
correct number of quark species, we use the fourth-root 
trick, which automatically performs the transition from 
four tastes to one taste per flavor for staggered fermion 
at all orders. We employ an interpolation operator with 
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isospin 7 = and = 0+ at the source and sink, 



where g is the index of the taste rephca, is the number 
of the taste rephcas, a is the color index, and we omit 
Dirac-spinor index. The time shce correlator C{t) for a 
meson in the momentum p state can be evaluated by 



onal sector 



+ (u^,(x,t)4(x,i)d^(0,0)d»(0,0)) 



(33) 



where we consider u, d quarks are degenerate, and let 0, x 
be the spatial points of the a state at source, sink, respec- 
tively. After performing Wick contractions of fcrmion 
fields, and summing over the taste index, for the light 
quark Dirac operator M , we obtain 

C(p, t)^-^Yj (Tril/^^ (x, t; x, t)TrM-\0, 0; 0, 0)) 

X 

+ E(-)"^'''" (Tr[Af-i(x, t; 0, 0)Af-i'x, t; 0, 0)]) , (34) 



where Tr is the trace over color index, and x = (x, t) is 
the lattice position. The first and second terms are the 
point-to-point quark-line disconnected contribution and 
connected contribution, respectively. 

For staggered quarks, the meson propagators have the 
generic single-particle form. 



Cit) 



i 



(35) 

where the oscillating terms correspond to a particle with 
opposite parity. For the a meson correlator, we consider 
only one mass with each parity in Eq. psp . namely, in 
our concrete calculation, our operator is the state with 
spin-taste assignment / (K) / and its oscillating term with 
7o75 ^ 7o75 [Ij, m, H^. Therefore, the a correlator was 
then fit to the following physical model, 

a(t) = 6,e-™'* + b„A-Ye~^'-Ai + {t^Nt-t), (36) 



Off- 



To avoid the complicated Fierz rearrangement of the 
quark lines, we choose the creation operators at the time 
slices which are different by one lattice time spacing as it 
is suggested in Ref. [s^, namely, we select ti = 0,t2 = 1, 
and = t for tttt — > a three-point function, and choose 
ti = 0, <2 = t, and ts = t + 1 for a — >■ tttt three-point 
function. 



(3- 





FIG. 2: Diagrams contributing to tttv — >■ a and a — >■ tttt three- 
point functions. Short bars stand for wall sources. Open cir- 
cles are sinks for local pion operator, (a) Quark contractions 
of TTTT — >■ a and (b) Quark contractions of a — ^ tttt. 

The quark line diagrams contributing to the a — > tttt 
or TTTT — >■ a three-point function denoted in Eq. (j24p are 
plotted in Fig. ^a.) and Fig. Eljb), respectively. The 
TTTT — a three-point function can be easily evaluated by 
constructing the corresponding three-point amplitudes 
for arbitrary values of the time slice using only two 
wall sources placed at the fixed time slices ti and ^2- 
However, the calculation of cr —>■ tttt three-point function 
is almost the same difficult as that of the rectangular di- 
agram for TTTT four-point correlator function, because it 
requires another quark propagator connecting the time 
slices t2 and ^3. For the nonzero momentum, we used 
an up quark source with 1, and an anti-up quark source 
with e'^P'^ on each site for pion creation operator. The 
a — >■ TTTT and tttt — >■ a three-point functions are schemat- 
ically shown in Fig. [TJ and we can also express them in 
terms of the quark propagators G, namely. 

Re J2 e'^"' (Tr[G,,(x3,t3)G^X3,i3)Gl^(xi,ti)]), 

Go-^7r7r(P, ^2, il) ~ 



where the 6,,^ and 6^ are two overlap factors. In Figure |3] Re ^ e*P ''^/Tr[Gtj (x2, ^2)^^^ (x2, ^2)^^^ (xi, ti)] \ ,(37) 



we can clearly note this oscillating term. 

We should understand that, for the staggered Kogut- 
Susskind quark action, our a interpolator couples to var- 
ious tastes as we examined in our previous study for the 
scalar ao and a mesons [H, [s^ , where we investigate 
the two-pseudoscalar intermediates states (namely, bub- 
ble contribution) . Therefore, we should remove a number 
of the unwanted two-pseudoscalar intermediates states 
with different tastes and slightly different energies from 
the (T correlator. 



D. Extraction of energies 

Through calculating the matrix of correlation function 
denoted in Eq. (p4)) . we can separate the ground state 
from first excited state in a clean way. It is very impor- 
tant to map out "avoided level crossings" between the a 
resonance and its decay products in a finite box volume. 
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because the first excited state is potentially close to the 
ground state. This makes the extraction of the ground 
state energy unfeasible if we only utilize a single exponen- 
tial fit ansatz. Since we can not predict a priori whether 
our energy eigenvalues are near to the resonance region 
or not, we find it safe to adopt the correlation matrix 
to analyze our lattice simulation data. To extract the 
lowest two energy eigcnstates, we follow the variational 
method [s^] and construct a ratio of correlation function 
matrices as 



Mit,tR)^Cit) C-\tR), 



(38) 



with some reference time slice 
to be large enough such that 
relation matrix M{t,t[i) from 
neglected, and the lowest two 
correlation function. The two 
be extracted by a proper fit to 
(n = 1,2) of matrix M{t,tfi) 
on the staggered fermion, and 
^n{t,tR) [n = 1,2) behaves as 

Xnit^tR) = An cosh 

(-l)*B„cosh 



tR [42[, which is assumed 
the contributions to cor- 
the excited states can be 
eigenstates dominate the 
lowest energy levels can 
two eigenvalues A„(t,tfl) 
Because here we work 
we can easily prove that 



(39) 



which explicitly contains an oscillating term, for a large 
t, which mean <^tR < i ^ T/2 to suppress the excited 
states and the unwanted thermal contributions. Here 
non-degenerate eigenvalues Xi{t,tR) > X2{t,tFi) are as- 
sumed. In practice, the oscillating term in Xi{t,tR) is 
not appreciable, we can also adopt following simple fit- 
ting model, 



Xi{t, tR) — A cosh 



-E t 



(40) 



and the difference between the fitting models of Eq. ([5^ 
and Eq. (|40]) is smah. 



III. LATTICE CALCULATION 
A. Simulation parameters 

We used the MILC lattices with 2 + 1 dynamical flavors 
of the Asqtad-improved staggered dynamical fermions, 
the detailed description of the simulation parameters can 
be found in Refs. [13, H^. We analyzed the tttt four- 
point functions on the 0.15 fm MILC lattice ensemble 
of 360 16^ X 48 gauge configurations with bare quark 
masses aniud = 0.0097 and arus = 0.0484 and bare 
gauge coupling lO/g^ = 6.572, which has a physical vol- 
ume approximately 2.5 fm. The inverse lattice spacing 
= 1.358t?^ GeV [Hill. The mass of the dynam- 
ical strange quark is near to its physical value, and the 
masses of the u and d quarks are degenerate. Periodic 
boundary condition is applied to three spatial directions 
and temporal direction. 



B. Sources 

To calculate the tttt correlation functions, we use the 
standard conjugate gradient method (CG) to obtain the 
required matrix element of the inverse fermion matrix. 
The calculation of the correlation function for the rect- 
angular diagrams naturally requires us to compute the 
propagators on all the time slices ts = 0, • • • , T — 1 of 
both source and sink, which requires the calculation of 
48 separate propagators in our lattice simulation. After 
averaging the correlator over all 48 possible values, the 
statistics are greatly improved since we can put the pion 
source at all possible time slices, namely, the correlator 
Cii(t) is then calculated through 



Cn(t) - ((7r^)(t)(w)t(0)^ 

i^((w)(i + t,)(w)t(t,,)), (41) 



The best-effort to generate propagators on all the time 
slices allows us to obtain the correlators with high pre- 
cision, which is very important to extract the desired 
scattering phase shifts reliably. 

For each time slice, six fermion matrix inversions are 
required corresponding to the possible 3 color choices for 
two pion sources, and each inversion takes about one 
thousand iterations during the CG calculation. There- 
fore, all together we carry out 288 inversions for each 
valence quark mass on a given configuration. As shown 
follow, this large number of inversions, performed on 360 
configurations, provides the substantial statistics needed 
to resolve the real part of the / = amplitude with reli- 
able accuracy. 

In the calculation of the off-diagonal correlator, C2i(t), 
the quark line contractions results in a three-point dia- 
gram. Since in this three-point diagram two pion fields 
are located at the source time slice ts, + 1, respectively. 
We calculate off-diagonal correlator C21 {i) through 

C2i(t) = (a(t)(7r7r)t(0)) = i ^ (a(t + ts)(^^)nts)) , 

(42) 

where, again, we sum the correlator over all time slices tg 
and average it. As for the second off-diagonal correlator 
Ci2(t), two pion fields are placed at the sink time slices 
ts + t and t + ts + 1, respectively, which make the com- 
putation of Ci2(t) difficult. However, using the relation 
Ci2{t) ~ C2i{t), we can obtain the off-diagonal matrix 
element C12 for free [Zsl . l46j . The a tttt component 
agrees with tttt — > a within the error, but the statistical 
errors of the matrix element C12 should be larger than 
that of matrix element C21 for a large t. Therefore, in 
the following analysis we substitute matrix element C12 
by the complex conjugate of matrix element C21, which 
is not only to save about 20% computation time, but also 
significantly to reduce statistical errors. 

For the a correlator, C22(t), we have measured the 
point-to-point correlator with high precision in our pre- 
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vious work for the this lattice ensemble, therefore, we 
can just use the available propagators to construct the 
CT-correlator 



(43) 



where, also, we sum the correlator over all time slices 
ts and average it. We must stress that we use Z2 noisy 
estimators based on random color fields to measure the 
disconnected correlator [Z^] . To obtain the signals prop- 
erly, we carry out numerical simulations with a high pre- 
cision and careful analyses. In our previous work |47|, we 
present the detailed procedures for the implementation of 
the Z2 method. Using the method discussed in Ref. (isj . 
we determine that 1000 noise Z2 sources are sufficiently 
reliable. 

One thing we must keep in mind that, in the calcu- 
lation of the correlator ((7r7r)(7r7r)^^), we make our best- 
efforts to reliably measure these rectangular and vacuum 
pieces, since the other two pieces are relatively easy to 
evaluate. We find that the rectangular and vacuum dia- 
gram make a major role in this correlator. Therefore, we 
get it properly for the tttt sector. 

In this work, we also measure two-point correlator for 
pion, namely, 



G.(t;p) = 



|7rt(O,t)7r(O,t,)|0), 
|7rt(p,t)^(p,i,)|0). 



(44) 



where the GTr{t;0) and GTr{t;p) is correlation function 
for the pion with the momentum and p, respectively. 



IV. SIMULATION RESULTS 
A. Time correlation function 

In Fig. [3] the individual ratios, which are defined in 
Eq. (|32|) corresponding to the diagrams in Fig. [l] 
{X = D,C, R and V) ^ are plotted as the functions of 
t. It is quite noisy for the disconnected diagram {V), 
but we can still get a good signal up to time separa- 
tion t = 14. The calculations of the amplitudes for the 
rectangular and vacuum diagrams stand for our princi- 
pal works. Clear signals observed up to t = 20 for the 
rectangular amplitude and up to t = 14 for the vacuum 
amplitude demonstrate that the method of wall source 
without gauge fixing used here is practically applicable. 

The values of the direct amplitude R^ is quite close 
to unity, indicating that the interaction in this channel 
is very weak. The crossed amplitude, on the other hand, 
increases linearly, which implies a repulsion in this chan- 
nel. After an initial increase up to t ~ 4, the rectangular 
amplitude exhibits a roughly linear decrease up to t ~ 15, 



^ We can easily verify that when t <C T/2, wc can approximately 
estimate the energy shift 5E from the ratio R-'^ . 
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FIG. 3: (color online). Individual amplitude ratios R''^ (t) for 
TTTT four-point function calculated by wall source as functions 
of t. Direct diagram shifted by 0.8 (diamonds), crossed dia- 
gram (octagons) and rectangular (squares) diagrams. 



and loss of signals after that, which suggests an attrac- 
tive force between two pions. These features are what we 
eagerly expected from the theoretical predictions [36l.l49j. 
We can observe that the crossed and rectangular ampli- 
tudes have the same value at i = 0, and the close values 
for small t. Because our analytical expressions in Eq. ((29)) 
for the two amplitudes coincide at t = 0, they should be- 
have similarly until the asymptotic tttt state is reached. 

The vacuum amplitude is negligibly small up to t ^ 
8 — 14, and loss of signals after that. This characteristic 
is in well accordance with the Okubo-Zweig-Iizuka (OZI) 
rule and xPT in leading order, which expects the disap- 
pearing of the vacuum amplitude [s^l • Since its correla- 
tion functions reflect the quantum fluctuations of QCD 
vacuum, the errors for disconnected amplitudes should 
be approximately independent of time separation t, and 
hence increases exponentially like e^™"*. Therefore, the 
signals present in this correlation functions are quickly 
buried in the noise and it is very difficult to obtain proper 
information from large time separation. 

In Fig. 131 we show the real parts of the diagonal com- 
ponents (tttt — >■ TTTT and a ^ a) and the real part of the 
off-diagonal component tttt — > cj of the correlation func- 
tion C{t) in Eq. ^I^. Since C{t) is a Hcrmitian matrix, 
hence, in the following analysis we substitute a — > tttt by 



9 



H H 



b e< 



~i — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — r 



G[ti7t— >a] 
□ G[ti7t:— >7nT] 




— I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — r- 
2 4 6 8 10 12 14 16 18 20 22 24 



<| 



10" 



10" 



10" 



O 

□o 



o Xi(t,tR) 

□ Xg(t,tR) 




10 

t 



15 



20 



FIG. 4: (color online). Real part of the diagonal compo- 
nents (tttt — >■ TTTT and a a) and the real part of the off- 
diagonal component nn a of the time correlation function 
C{t). Occasional points with negative central values for the 
off-diagonal component tttt — >■ a and the diagonal component 
CT — >■ (T are not displayed. 

TTTT — > (T to reduce errors. 

As wc pointed out in Sec. Ill C 21 there exists bubble 
contribution in sigma correlator. In this work we will 
calculate the scattering phase shift with the bubble con- 
tribution removed from corresponding sigma correlator. 
In Ref. [1^ H^, we presented the bubble term Bi^{t), 
which was parameterized by three low-energy couplings 
5A^ and 5v- In our concrete calculation, they were 
fixed to the values of our previous determinations f23l . l40| . 
The taste multiplct masses in bubble terms were fixed as 
listed in Table 2 in Ref. [2^ . After deducting the bubble 
term, the remaining sigma correlator now contains the 
clean information for sigma meson. 

Wc calculate the two eigenvalues A„(t,tfl) (n — 1,2) 
for the matrix M{t,tfj) in Eq. ([55)1 with the reference 
time tji = 6. In Fig. [S] we plot our lattice results for 
^n{t,tfi)(n = 1,2) in a logarithmic scale for the moving 
frame, as a function of time t together with a correlated 
fit to the asymptotic form given in Eq. pop . From these 
fits we then extract the energies that will be used to 
determine the scattering phase. 

To extract the energies reliably, we must take two ma- 
jor sources of systematic error into consideration. One 
arises from the excited states which affect the correla- 
tor in low time slice region. The other one stems from 
the thermal contributions which distort the correlator in 
high time slice region [s^]. By denoting a fitting range 
[tmin,imax] and varying the values of the iniin and imax, 
we can control these systematic errors. In our concrete 
fitting, we take tmin to be -|- 1 and increase the refer- 
ence time slice tu to suppress the excited state contam- 
inations. Moreover, we select tmax to be sufficiently far 
away from the time slice t = T/2 to avoid the unwanted 



FIG. 5: (color online). The eigenvalues Xi{t,tii) and \2{t,tR) 
as the function of t are shown. The solid lines are correlated 
fits to Eq. (|39p . Occasional points with negative central values 
for the correlator \2{t,tfi) axe not plotted. 



thermal contributions. The corresponding fitting param- 
eters in, tmin and tinax uscd in this work are tabulated 
in Table IH The corresponding fittings result in the rea- 
sonable values of x^/dof. The x^/dof together with the 
fit results for En (n = 1, 2) are also listed in Table ID 



TABLE I: The values of the energy eigenvalues for the ground 
state (n = 1) and the first excited state (n = 2). In the table 
we list the reference time Ir, the lower and upper bound of 
the fitting range, tmin and tmax, the number of degrees of 
freedom (dof) for the fit quality x^/dof and the fit results for 
the energy eigenvalues En (n = 1, 2) in lattice units. 



n 


tR 


^min 


imax 


aEn 


xVdof 


1 


6 


7 


12 


0.6767(20) 


0.922/2 


2 


5 


6 


12 


0.8086(75) 


0.158/3 



The energy of the free pions Ei is calculated from the 
mass and energy E-,^ obtained by a single exponential 
fit to G^(t;0) and G^(t;p) in Eq. (gl]), as Ei=mT, + E^. 
These results are summarized in the upper part of Ta- 
ble ini We observe that Ei < Ei < E2, which mean 
that the phase shift for Xi{t,tji) and X2{t,tji) is posi- 
tive and negative, respectively. This clearly indicates the 
existence of a resonance in between. 

In Table Hill we show the mass m and the energy E of 
the pion and a meson with momentum p = (27r/L)e3, 
calculated from the corresponding time correlation func- 
tions. The a meson with the zero momentum can not 
decay energetically, therefore, the mass and energy are 
extracted in a usual way [23j . 
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TABLE II: The energy of tttt system En and the scatter- 
ing phase shift S. Ei is the energy of the free pions. E„ is 
obtained from the fitting. The invariant mass ^/s, the the mo- 
mentum k and the phase shifts S calculated with the energy- 
momentum relations (|47fl in the continuum are referred to as 
Cont, and those obtained with the relations (|48p on the lat- 
tice are referred to as Lat. The momentum fco is defined by 
fco = s/4 — m^. All values with the mass dimension are in 
lattice units. 





n = 1 




n = 2 




En 


0.7085(6) 








En 


0.6767(20) 




0.8086(38) 






Cont 


Lat 


Cont 


Lat 




0.5511(25) 


0.5613(25) 


0.7068(43) 


0.7179(44) 




0.0155(7) 


0.0185(8) 


0.0644(15) 


0.0699(17) 


kl 




0.0182(7) 




0.0684(16) 


tan (5 


0.380(17) 


0.257(14) 


-1.261(44) 


-1.509(55) 


sin^ 5 


0.126(9) 


0.0621(64) 


0.614(16) 


0.695(15) 



Therefore, the definitions of y/s and k contain the similar 
discretization errors, which have been extensively studied 
by Rummukainen and Gottlieb [s^], and they suggest 
using the lattice modified relations, namely, 

cosh(-y/s) ~ cosh(_E7i/i?) — 2 sin^(p/2), 

2sin2(fc/2) = cosh - cosh(m^), (48) 

and hence the scattering phase shift is obtained by plug- 
ging the momentum k into the formula in Eq. (|2ip. 

To understand these discretization effects, in the 
present study we calculate the invariant mass -y/s and 
the momentum k from the energy momentum relations 
both in the continuum (|^7|) and on the lattice (1^5]) . and 
then estimate the corresponding phase shift by inserting 
the k into Eq. ((2T|) . We refer to the difference stemming 
from two choices as the discretization error, which should 
disappear in the continuum limit. The results for the in- 
variant mass -^/s (or Ecm), the momentum k and the 
phase shift 5 are tabulated in Table |IT1 



TABLE III: Mass m of the tt and a meson, and energy E 
of the TT and a meson with the momentum p = (27r/L)e3, 
extracted from the time correlation function. All values are 
in lattice units. 





TV 


a 


m 


0.2459(2) 


0.594(35) 


E 


0.4626(5) 


0.714(22) 



B. Lattice discretization effects 

In the continuum limit we have a dispersion relation 

E = \/ p'^ + rn?, 



(45) 



for the single particle, which is particularly relevant for 
the finite-size methods used here |34ll. On the lattice, this 
relation should be modified to 1341 



cosh{E) = 2 • sin^(p/2) + cosh(m). 



(46) 



We should consider the important discretization error 
in Rummukainen- Gottlieb formula ()2ip . It arises from 
the Lorentz transformation from the moving frame to 
the center-of-mass frame using Lorentz symmetry in the 
continuum limit. In the Lorentz transformation we use 
these relations. 



(47) 



for the invariant mass y/s, the energy in moving frame 
Emf and the momentum k. However, on the lattice, the 
discretization effects explicitly break Lorentz symmetry 
and Eq. (|47p is only valid up to the discretization errors. 



C. Extraction of resonance parameters 

From Table HIl we note that the considerable differ- 
ences which depend on the chosen energy-momentum re- 
lations are obviously observed in ^/s and k. Moreover, 
the thence difference for phase shift 5 is significantly 
larger than the corresponding statistical errors. These 
are also shown in Fig. |6l where the phase shift sin^ 5, 
which is proportional to the scattering cross section of 
TTTT system, is drawn. In Table |ll] we see that the sign of 
the phase shift 5 at y/s < nia- {ama = 0.594(33)) is pos- 
itive, which indicates an attractive interaction, and that 
at -^/s > rria is negative, which suggests a repulsive inter- 
action. These features are what we expected. It confirms 
the presence of a resonance at a mass around the <t mass 

In principle, it is straightforward to extract the a me- 
son decay width by fitting the phase shift data with the 
RBWF. However, the quark mass we studied here is sig- 
nificantly larger than its nature value, therefore a long 
chiral extrapolation is required. Since the kinematic fac- 
tors in the decay width depend explicitly on the quark 
mass [2^ . and we just made a lattice simulation on one 
set of the quark mass, we have to avert this direct mea- 
surement of the decay width and adopt an alternative ap- 
proach. As we discussed in Section fll Al we parametrize 
the resonant characteristic of the s-wave phase shift i5o 
in terms of the effective cr — >■ tttt coupling constant Qa-mr, 
namely. 



tan i5o 



2 



Sir ^s{Ml-sy 



where Mr is the resonance mass and ga 
through the lagrangian 

LcS ^ ga-UTT CTTTTT, 



(49) 



is denoted 



(50) 
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as in the continuum theory. According to the discus- 
sions in Ref. [22|, we can reasonably assume that the 
couphng constant ga-Tnr vary rather slowly as the quark 
mass changes. Thus, the equation (|49p allows us to solve 
for two unknown parameters, namely, the coupling con- 
stant gaTTir, and the a meson mass m„. 

The invariant mass ^/s and the momentum k appear- 
ing in the continuum effective theory in Eq. ()49|) satisfy 
the energy- momentum relations ()47|) . The thence dis- 
cretization error may arise in the choice of ^/s and k in 
the application of Eq. (|49)) to the extracted phase shift. 
Fortunately, our lattice results show that this does not 
cause a serious problem numerically. In Table we also 
present the momentum evaluated by /cq = s/4 — 
for a given y^, and m^. We can notice that the difference 
between k and ko is not remarkable. Therefore, the error 
stemmed from the different definitions of the scattering 
momentum can be reasonably ignored. In practice, we 
employ the scattering momentum fcg upon the applica- 
tion of Eq. (PI). 

The lattice simulation results of the coupling constant 
ga-TTTT and the resonance mass Mr solved by Eq. are 



g,^^ = 3.22(52) GeV, 
Mr = 0.660(31), 
Majm, = 1.112(85), 



(51) 



where we utilize the energy-momentum relations (j47p in 
the continuum, and the a meson mass is extracted 
from the corresponding correlation function in our previ- 
ous study. If we employ the energy momentum relations 
()48p on the lattice, we obtain the corresponding simula- 
tion results as 



Mr 

MR/m„ 



2.69(44) GcV, 
0.691(37) , 
1.163(93) . 



(52) 



The value of coupling constant gaTm is in reason- 
able agreement with ga-wTr = 2.47(45) GeV obtained in 
Ref. [Hil, g„„^ = 2.97(4) GeV obtained in Ref. ^ and 
g„^^ = 2.86 GeV given in Ref. [22|. 

In Fig. [6l we display the curves for sin^ So obtained by 
Eq. (j49p with the coupling constant go-iTTr and the reso- 
nance mass Mr given in Eq. ((5T|) and Eq. (|52p . respec- 
tively. The position at sin^ (5o = 1 which stand for the 
resonance mass AIr is also displayed in Fig. [5] for two 
cases (black cross and red plus for the continuum and 
lattice case, respectively) . For visualized comparison, we 
also draw the sigma mass (fancy cyan plus). We can 
note that Mr is in reasonable agreement with the sigma 
mass TOct- 

Assurning that the dependence of gaTnr on quark mass 
is small [22| , we can roughly estimate the a meson decay 
width at the physical quark mass as 
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(53) 



FIG. 6: (color online). The scattering phase shift sin^ 5, po- 
sitions of a mass rricr and resonance mass Mr. Cont refer to 
the results achieved with the energy momentum relations in 
the continuum limit (|47p and Lat to those with the relations 
on the lattice (|48p . The two lines are obtained by Eq. (|49p 
with parameters ga-mr and Mr given in Eq. (|5ip and Eq. (|52l) . 
respectively. The abscissa is in lattice units. 



where m^^^ — 513 ± 32 MeV is the estimated physical a 
meson mass taken from PDG [l[ , and momentum k^^'^ is 
calculated by 

(^phy)2 ^ (^.^ ) _ (^phy)2 ^ 

where is physical pion mass {m}^^ w 140 MeV) [l[. 
This produces 

FP'^y = (337 ± 82) MeV (54) 

where we use the data given in Eq. (j5ip . and 

pphy = (236 ± 49) MeV (55) 

where we employ the data given in Eq. (|52p . We can ob- 
serve that the difference stemming from our two chosen 
energy-momentum relations is comparable with the sta- 
tistical error. Although our preliminary estimates for the 
a — !■ TTTT decay width in this work is not within the PDG 
estimated resuh P = 600 - lOOOMeV [l|, this is stiU an 
inspiring result, considering that we make a big assump- 
tion that the coupling constant does not depend on the 
quark mass, and we perform a long chiral extrapolation, 
etc. 

In the present study, we make an extensive use of the 
relativistic Breit Wigner formalism, and the sigma me- 
son is a very wide object and the RBWF approxima- 
tion holds perfectly for relatively narrower objects. As 
it discussed in Ref. [53[, we should adopt a much more 
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modcl-indcpendcnt approach to the extraction of the fi- 
nite volume hmit. In our future tasks, we must address 
the phenomenological treatment. 

Bearing in mind that this work is an exploratory study, 
the main purpose is to present some conceptual issues of 
determining o resonance parameters directly from lattice 
QCD calculations, wc think that it is enough justified all 
these above assumptions and simplifications. 

V. CONCLUSIONS AND OUTLOOKS 

In the present work, we have carried out a direct lattice 
QCD calculation of the s-wave tttt scattering phase shift 
for the isospin / = channel near the cr-meson resonance 
region at total non-zero momentum in the moving frame, 
where the rectangular graph and vacuum diagram play 
a key role, for the MILC "medium" coarse (a = 0.15 
fm) lattice ensemble in the presence of 2 -I- 1 flavors of 
the Asqtad improved staggered dynamical sea quarks. 
We employed the same technique in Ref. [13, , which 
we use the moving pion wall source operators without 
gauge fixing for the / = channel to obtain the reliable 
precision. We calculated all the four diagrams which are 
categorized in Ref. [13, [13 , and observed a clear signal 
of the attraction for / = channel. 

We have demonstrated that the calculation of the s- 
wave scattering phase shift for the / = tttt system and 
then the estimation of the decay width of a meson are 
practicable on the lattice QCD with our current limited 
computing resources. The phase shift data clearly shows 
the presence of a resonance at a mass around the a meson 
mass obtained in our previous study [23i] . This resonance 
can be reasonably identified with the a meson. Moreover, 
we extracted the a meson decay width from the phase 
shift data and showed that it is fairly compared with the 
estimated a meson decay width from PDG [l|. 

However, we realize some important issues which 
should be cleared in the future works. One is to reduce 
the discretization errors, which wc show in the previ- 
ous section are significantly larger than the correspond- 
ing statistical errors. A simple way to solve this problem 
is to use a lattice configuration closer to the continuum 
limit. Moreover, a full analysis to determine the size de- 
pendence of the phase shift employing a set of lattice 
sizes is also highly desired. However, all of these issues 



are beyond the scope of this paper since this will require 
a substantial amount of computing resources. We leave 
these challenging tasks in our future study. 

We adopted the effective range formula, which allows 
us to use the effective a — >■ tttt coupling constant ^o-tttt to 
extrapolate from our lattice simulation point m-^/ma w 
0.414 to the physical point mT^jma sa 0.273, assuming 
that ga-K-K does not depend on the quark mass. This is 
just a rough estimation, thus a more direct evaluation 
of the decay width is highly desirable. As we pointed 
above, the decay width can be estimated directly from 
the energy dependence of the phase shift data by fitting 
the BWRF if we make the lattice simulations near to 
the physical quark mass and we have simulation data 
which have several energy near the resonance mass as it 
was done for the calculation of the p meson decay width 
in Ref. [s^]- We will continue to look for the possible 
computational resources to accomplish this task. 

When our preliminary lattice results obtained are com- 
pared with the corresponding PDG quantities, it is clear 
that the lattice QCD simulations can not yet match the 
accuracy of PDG, and even can not be considered to 
be "physical" one. Although a precise determination of 
the cr resonance parameters on the lattice is absolutely 
a big challenge, our preliminary work reported here can 
serve as stepping out the first step in an attempt to un- 
derstand the strong a decay from direct lattice QCD in 
a conceptually clean way. 
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